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ABSTRACT 

The evolution of dust at redshifts z > 9, and consequently the dust properties, differs greatly 
from that in the local universe. In contrast to the local universe, core collapse supernovae 
(CCSNe) are the only source of thermally-condensed dust. Because of the low initial dust-to- 
gas mass ratio, grain destruction rates are low, so that CCSNe are net producers of interstellar 
dust. Galaxies with large initial gas mass or high mass infall rate will therefore have a more 
rapid net rate of dust production comported to galaxies with lower gas mass, even at the same 
star formation rate. The dust composition is dominated by silicates, which exhibit a strong rise 
in the UV opacity near the Lyman break. This “silicate-UV break” may be confused with the 
Lyman break, resulting in a misidentification of a galaxies’ photometric redshift. In this paper we 
demonstrate these effects by analyzing the spectral energy distribution (SED) of MACS1149-JD, 
a lensed galaxy at z = 9.6. A potential 2mm counterpart of MACS1149-JD has been identified 
with GISMO. While additional observations are required to corroborate this identification, we 
use this possible association to illustrate the physical processes and the observational effects of 
dust in the very high redshift universe. 

Subject headings: galaxies: high-redshift - galaxies: evolution - galaxies: individual (MACS1149- 
JD) - Interstellar medium (ISM), nebulae: dust, extinction - physical data and processes: nuclear 
reactions, nucleosynthesis, abundances 


1. INTRODUCTION 

Galaxies are the principle tracers of the global evolution of the star formation history, the metal en- 
richment, and the formation of dust in the universe. Their detection at very high redshifts offers a unique 
glimpse of these physical processes during the earliest epochs of galaxy evolution. Dust has been detected 
in several high-redshift galaxies (Gall et al. 2011b), the current record being held by the dusty quasar QSO 
J1120+0641 located at redshift z = 7.12 (Venemans et al. 2012). Dust attenuates the intrinsic spectrum 
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of the stellar population, affecting the various physical quantities derived from UV-optical observations. 
Detection of dust through its extinction and emission is therefore essential in determining the properties of 
galaxies in the high-redshift universe. 

A gravitationally magnified galaxy, designated MACS1149-JD was recently discovered at a photometric 
redshift of z = 9.6±0.2 (Zheng et al. 2012). The massive lensing cluster (MACS J1149. 6+2223), located 
at redshift 0.544, provides a magnification factor of /i = 14.5±1 4 'q. Assuming no extinction by dust, and 
correcting for magnification, the star formation rate (SFR), the total stellar mass, and luminosity derived 
from the observed UV-optical (UVO) flux are about 1.2 /z^ 1 M 0 yr -1 , 1.5 x 10 8 /z)^ 1 M 0 , and 2 x 10 10 /z^ 1 L & , 
respectively, where /Z15 is the lensing magnification normalized to the nominal value of 15 (Zheng et al. 2012). 
However, uncertainties in fitting the spectral energy distribution (SED) of this galaxy allow the possibility 
of a certain amount of extinction (Zheng et al. 2012). 

The Goddard-IRAM Superconducting 2 Millimeter Observer (GISMO) camera at the IRAM 30m tele- 
scope on Pico Veleta (Staguhn et al. 2013), Spain was used to probe whether MACS1149-JD has a rest frame 
far-infrared counterpart. The observations detected a 2 mm source with an intensity of 400 ± 98 /zJy (see 
Figure 1). The observed offset of 8" (less than half the GISMO beam) is entirely consistent with this source 
being at the HST detected position of MACS1149-JD. The rest wavelength is ~ 190 /zm, close to the peak of 
the dust emission in star forming galaxies. These observations strongly suggest that dust is present in this 
galaxy, however, further confirmations are needed to absolutely corroborate this association. Observations 
using the SHARC-2 submillimeter bolometer camera at the Caltech submillimeter observatory yielded a 
350 /zm upper limit of 8.6 ± 5.6 mJy within an area consistent with the nominal position of MACS1149-JD. 

For sake of the current analysis we will assume the validity of this association, adopt the observed 
2 mm flux as originating from MACS1149-JD, and demonstrate the unique properties of any dust that 
may be forming in the high redshift universe. All fluxes and results of this paper are presented in the 
galaxy’s rest frame with the adopted cosmological parameters Ho = 67 km s -1 Mpc -1 , Om = 0.315, 
and Ua = 0.685 (Planck Collaboration et al. 2013), and a lensing magnification factor of 15. With these 
cosmological parameters the age of the universe is ~ 500 Myr at z = 9.6. 


2. MODELING THE SED OF MACS1149-JD 

In the absence of any starlight, dust would be in thermal equilibrium with the cosmic microwave back- 
ground (CMB). This effect can be important at high redshifts where starlight and the CMB can have 
comparable energy densities (da Cunha et al. 2013). The GISMO flux is the differential 2 mm flux with 
respect to the CMB emission. Therefore, the observed 2 mm flux consists of only reradiated starlight. We 
characterize the IR spectrum by T d , the temperature that would be attained by the dust in the absence 
of the CMB, and by a mass absorption coefficient k(A), that has a X^ 13 (/3 = 1.5) dependence at far-IR 
wavelengths and a value of 1.5 cm 2 g -1 at A = 850 /zm (Kovacs et al. 2010). The specific IR luminosity in 
the galaxy’s rest frame is then given by: 

L„(A) = 4M d «(A) [7rR,(A,T e// )-7rR,(A,T cmb )] (1) 

where Teff = (T d 4+0 + T^) 1 ^, M d is the dust mass, and T cmb = (1 + z) 2.73 ~ 29 K is the CMB 
temperature at z = 9.6. L„( A) is a function of the dust temperature and mass. With L„( A = 190/zm) 
normalized to the observed GISMO flux, the dust temperature T d and dust mass are inversely related. 


The presence of dust significantly affects the properties of AIACS1149-JD derived under the assumption 
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that it is dust free, and constrains viable star formation scenarios for this galaxy. The panels of Figure 2 
depict the analysis used to constrain the star formation history (SFH), the amount of dust obscuration and 
resulting infrared emission, and the mass of radiating dust, based on the available UVO data and adopted 
far-IR flux. Two star formation histories, characterized by a delayed exponential star formation rate (SFR), 
ip{t) = ipo ( t/t exp ) exp(— t/t exp +l), were used in the analysis with values of t exp =100 and 300 Myr (hereafter 
models t exp 100 and t exp 300, respectively, see Figure 2a). 

Intrinsic stellar spectra were calculated using the stellar population synthesis code PEGASE (Fioc & 
Rocca-Volmerange 1997) with low-metallicity (LowZ) and Kroupa (Kr) stellar initial mass functions (IMFs) 
(Marks et al. 2012; Kroupa 2001). The LowZ IMF is more heavily weighted towards the formation of massive 
stars, and its stellar population produces more UV radiation compared to its Kroupa counterpart. Its choice 
is motivated by the theoretical argument that inefficient cooling inhibits the fragmentation of protostellar 
clouds into smaller units (Bromm & Larson 2004), and by observed variations of the IMF with stellar 
metallicity (Marks et al. 2012; Kroupa et al. 2013). Figure 2b shows the spectra for ip o = 1 Mq yr” 1 . 

For each star formation (SF) scenario, characterized by t exp , ipo, and a stellar IMF, we calculated the 
amount of dust absorption needed to reproduce the observed UVO fluxes. The ratio between the observed 
and intrinsic flux is the escape probability at that given wavelength, which for a homogeneous sphere is given 
by (Lucy et al. 1991; Cox & Mathews 1969; Osterbrock & Ferland 2006): 

P esc (r,w) = P e5C (r, 0)[l-cu + o;P esc (T, 0)] (2) 

P esc (r,0) = A [1 - + (- + ^L) exp(-2r)] 

where t(A) = (3/4)Md K(X)/nR 2 is the optical depth, and w( A) is the dust albedo. The second expression 
corresponds to the case when w = 0. The albedo of small dust grains (a < 0.01 /xir i) is < 0.2 and P esc (r, 0) ~ 
0.1 for most scenarios. With these parameters, scattering will increase the escape probability by at most 
20%, having only a minor effect on the conclusions of this paper. 

We assumed that the IR emission is only powered by stellar radiation, so that the dust giving rise to 
the IR emission is also responsible for the attenuation of the stellar emission. The reradiated IR luminosity, 
given by the integral of L„( A) in eq. (1), must then be equal to that the starlight energy absorbed by the 
dust: 

Labs = J Ll nt (X)[l-P esc (X)]du (3) 

where L(, nt is the intrinsic stellar luminosity. Figure 2c illustrates this principle of energy balance for model 
t exp 300-Low Z for values of ipo = 20 and 40 Mq yr” 1 . The dust must correspondingly radiate at temperatures 
of 44 and 70 K, with dust masses of 1.2 x 10' and 4.0 x 10 6 Mq, respectively. 

Not all SF scenarios can fulfill the energy constraint. Figure 2d plots the IR luminosity as a function of 
Td , and the absorbed stellar luminosity as a function of ipo- Energy conservation requires Lir (red curve) 
to be equal to L a j, s (green and blue curves). The GISMO detection sets a lower limit on the dust luminosity 
regardless of the dust temperature. This constraint sets a corresponding lower limit on ipo for each SF 
scenario (see Table 1 for details). 

An additional constraint on viable SF scenarios is that the dust mass giving rise to the UVO attenuation 
and the IR emission be produced within the allotted time of < 500 Myr, taking possible grain destruction 
processes into account (Dwek et al. 2007). 
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3. THE EVOLUTION OF DUST AT HIGH REDSHIFTS 

Dust in galaxies is produced in the ejecta of explosive core-collapse supernovae (CCSNe) and in the 
quiescent winds of asymptotic giant branch (AGB) stars which are not massive enough to end their life 
as CCSNe. Refractory elements locked up in dust are returned to the gas phase by thermal and kinetic 
sputtering and evaporative grain-grain collisional processes in supernova blast waves (Jones 2004). Star 
formation removes both, gas and dust from the ISM. The equations governing the evolution of the dust and 
gas masses in the ISM of galaxies were presented in Dwek & Cherchneff (2011). At high redshifts, CCSNe 
dominate the rate of dust production at a rate given by dMd/dt = YdRsN , where Yd is the IMF-averaged 
dust yield in CCSNe (see Dwek et al. 2007, Table 2 for the maximum attainable dust yield as a function 
of stellar mass), and Rsn oc tpo is the IMF-dependent rate of CCSNe. CCSNe are also the main source of 
grain destruction in the ISM, at a rate given by Md/rd , where Td(t) is the dust lifetime, given by (Dwek & 
Scalo 1980): 

T _ Md{t) _ Mgjjf) 

d fn d {t)RsN(t) frig Rsn (t ) 

where ffid(t) is the time-dependent total mass of refractory elements locked in dust that is returned back to 
the gas phase by a single supernova remnant, fn g = frid(t)/Zd(t) is the effective mass of interstellar gas that 
is totally cleared of dust, where Z d (t) = Md/M g is the dust-to-gas mass ratio in the ISM, and M g is the 
total gas mass. The dust destruction rate is then given by: 

= fn d {t) RsN(t) = Z d {t)fn g Rswit) (5) 


The value of fn g is approximately time-independent and about equal to 10 3 M 0 (Dwek et al. 2007), and 
Zd « 0.007 in the solar neighborhood (Zubko et al. 2004). The total mass of dust destroyed by a single SNR 
is therefore ~ 7 M 0 , significantly larger than the total ~ 1 M 0 of condensible elements in their ejecta. In 
the local universe, CCSN are therefore net destroyers of interstellar dust. 

In a dust-free galaxy, CCSNe are obviously net sources of interstellar dust. There is therefore a limiting 
value of Zd ~ 1 x 10~ 4 for which the grain destruction rate by SNR is balanced by their formation rate, 
assuming a 10% condensation efficiency in the ejecta. The initial rate of dust enrichment in a galaxy depends 
therefore on: (1) the stellar IMF - a top-heavy IMF will have a higher IMF-averaged dust yield; and (2) the 
evolution of the gas mass in the galaxy, which in the instantaneous recycling approximation is given by: 


Mg(t) = Mg{ 0) 


(m e j) 

(m) 


ip{t') dt' 



(6) 


where M g ( 0) is the initial gas mass, (m e j) and (to) are, respectively, the IMF-averaged ejecta and initial 
stellar mass, and ( dM g /dt)i n f is the infall rate of the gas into the galaxy. 

Dust masses derived from evolutionary models depend on tpo and frig, whereas dust masses derived from 
fitting the UVO-IR SED depend only on ip o- Figures 2e-f depict these two distinctly derived dust masses as 
a function of the stellar mass derived from the corresponding population synthesis and dust evolution models 
at t = 500 Myr, the epoch of observations. Figure 2e represents the results for the LowZ IMF and an initial 
gas mass of 1 x 10 10 M 0 , and Figure 2f for the Kroupa IMF and an initial gas mass of 5 x 10 10 M 0 with no 
infall. The solid black lines show the relation between Md and stellar mass for the dust evolutionary models 
with selected lines labeled by the value of frig (in M 0 ). Tick marks along the lines indicate the values of ipo- 
The models were calculated assuming that all refractory elements condensed in the mass outflow from AGB 
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stars, but that only 10% did in the ejecta of CCSN. The thick dashed lines show the derived dust masses in 
the absence of any grain destruction. The rate of grain destruction depends on the dust-to-gas mass ratio 
in the ISM, and therefore on the choice of the initial mass of the gas reservoir. The line connecting the 
filled and open diamonds represents the relation between the dust mass derived from the SED fitting to the 
stellar mass derived from the population synthesis models. The diamonds are labeled by the value of i/'o- 
The hatched region gives the dust masses for all combinations of ipo and fn g considered in this paper. Viable 
SF scenarios are those for which dust masses derived from SED fitting are equal to those derived from the 
evolutionary models (hatched regions). The resulting constraints on the SF scenarios and dust and stellar 
masses are summarized in the 3rd and 4th column of Table 1. 

In general, SF scenarios characterized by a LowZ IMF produce the inferred dust mass at lower SFR 
and lower stellar masses. Scenarios that require low values of m g imply dust lifetimes that are significantly 
larger than those expected for a homogeneous ISM, implying a clumpy ISM (Dwek & Scalo 1979; Dwek et al. 
2007). 

Figure 3 depicts the evolution of the dust for different values of the SFR and the stellar IMF, highlighting 
the role of the gas mass in determining the dust production rate. The models assume no infall, so the gas 
mass is only determined by its initial value and net consumption rate by stars. Models with large initial gas 
masses produce dust at a higher rate than those with a lower initial gas mass. Eventually, the dust-to-gas 
mass ratio will reach the critical value of ~ 0.001 after which the rate of grain destruction by CCSN exceeds 
the rate of dust formation in their ejecta. Galaxies with large initial gas mass have therefore a faster rate 
of dust enrichment. However, they will also be rarer objects since lower mass galaxies are more readily 
assembled at z ~ 10 than the more massive ones (Somerville et al. 2008; Behroozi et al. 2013). 


4. THE ULTRAVIOLET SILICATE BREAK 

The galaxy’s intrinsic spectrum was forced to reproduce the observed UVO emission (see Figure 2c). It 
is therefore of interest to examine how the derived mass absorption coefficient of the dust compares with that 
of astrophysical dust grains. Figure 4 compares the wavelength dependence of k(A) with that of 0.01 /im radii 
silicate and amorphous carbon grains, with optical constants from Li & Draine (2001). The mass absorption 
coefficient shows a sharp rise at A < 0.15 /tm. The rise is similar to that of the silicate grains but at shorter 
wavelengths. A good match would require the galaxy to be at a somewhat lower redshift than that inferred 
from the UVO SED fitting. This suggests that young galaxies, in which silicates are the dominant dust 
component, will exhibit a sharp drop in their spectrum around 0.15 /im caused by this rising silicate UV 
opacity. This “UV silicate break” could be mistaken for an intrinsic Lyman break at A = 0.0912 /im, leading 
to an overestimate in the photometric redshift determination of galaxies. Such galaxies may also exhibit 
a break in their IR spectrum that is associated with the 9.7 and 18 /mi silicate absorption features. This 
“IR silicate break” can be used to photometrically select galaxies in the z ~ 1 — 2 range (Takagi & Pearson 
2005), or to estimate the photometric redshift of high-z galaxies. 


5. SUMMARY 

The evolution of dust at redshifts above ~ 10 differs fundamentally from that at lower redshifts: (1) 
at these redshifts, CCSNe are the dominant sources of interstellar dust. The enrichment of the ISM with 
carbon and carbon dust is delayed until the most efficient carbon producing AGB stars have evolved off the 
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main sequence. The composition of the ISM and its dust may therefore be significantly different from that 
in lower redshift galaxies (Dwek 1998). The redshift interval over which this effect can be observable may be 
quite narrow, depending on the stellar IMF and the mass range of carbon rich stars; (2) at low metallicities, 
hence low dust-to-gas mass ratios, CCSNe are net producers of interstellar dust. This is in sharp contrast to 
lower redshift galaxies, where CCSNe destroy more dust than they produce, which has led to the suggestions 
that most of the interstellar dust in galaxies must be reconstituted by cold accretion in the ISM (Dwek & 
Scalo 1980; Zhukovska et al. 2008; Draine 2009; Calura et al. 2010; Valiante et al. 2011; Gall et al. 2011a; 
Zhukovska 2014); (3) the efficient production of dust in the high-z universe is also greatly facilitated by 
the fact that at low-metallicities the stellar IMF is biased towards the formation of massive stars, further 
reducing the relative importance of lower mass AGB stars to the production of carbon in the early universe. 

Analysis of the SED of MACS1149-JD, including the 2 mm (190 pm restframe) flux from its GISMO 
counterpart, shows that dust formation proceeds at rapid rate at the earliest stages of galaxy evolution. 
Dust masses of ~ (2 — 7) x 10 7 Mq can be readily produced within 500 Myr of evolution with SFR as low as 
5 — 8 Mq yr -1 , without the need to resort to growth in the ISM. This dust exhibits a sharp rise in opacity 
at UV wavelengths around 0.15 pm, characteristic of very small silicate grains expected to form in the ejecta 
of CCSNe (Nozawa et al. 2010; Cherchneff & Dwek 2010). The sharp rise produces a ”UV silicate break” in 
the spectra of very young galaxies, which can be mistaken for a Lyman break, leading to an overestimate of 
their photometric-based redshifts. The existence of high-redshift, actively star-forming, galaxies with little 
dust, such as Himiko at £ = 6.62 (Ouchi et al. 2013), shows that the dust abundance in galaxies strongly 
depends on their star formation history, the efficiency of grain destruction processes, and the presence of gas 
infalls and dusty outflows. Our analysis shows that dust can form following the death of the first massive 
stars, but its survival and abundance will depend on the particular physical processes operating in each 
individual galaxy. Finally, independent confirmation of the association of the GISMO 2mm source with 
MACS1149-JD, will render this galaxy the youngest dust forming galaxy, with profound implications for the 
onset of radiative, thermal, and chemical processes in the early universe. 
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Fig. 1.— Left: The GISMO 2 mm signal-to-noise map of the MACS1149-JD field. The 4.1er source near the center 
of the field is is ~ 8" (or < 0.5 beam widths [FWHM] = 17. 5") from the Hubble location of MACS1149-JD (red 
cross), and therefore consistent with being the same source. Two other sources with S/N > 3 are present in the field. 
Right: HST WFC3 image of MACS1149-JD at 0.85, 1.05, and 1.60 /im (blue, green, red respectively). MACS1149- 
JD (arrowed) can be seen as the only extremely red source (a dropout at 0.85 and 1.05/rm). The field shown here 
is much smaller than the GISMO field. The blue circle indicates the 24. 7" FWHM of the smoothed GISMO beam 
centered on the nominal location of the GISMO source. 
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Fig. 2. — (a) The star formation rates as a function of cosmic age considered in this study have time scales of 
texp = 100 or 300 Myr. At the time when MACS1149-JD is observed (t = 500 Myr) the SFR has fallen to 0.092 or 
0.86 of the peak SFR, ipo] (b) the spectrum of the galaxy at t = 500 Myr for the different SFRs shown in (a) for the 
Low-metallicity and Kroupa IMF. (c) models for the intrinsic spectrum must satisfy the energy constraint, requiring 
the absorbed stellar luminosity to be equal to the reradiated IR Luminosity. The two examples depicted in the figure 
are characterized by values of ipo = 10 and 40 Mq yr _1 , and corresponding dust temperatures of T r j = 44 and 70 K; 
(d) depiction of SF scenarios that satisfy the energy constraint. The scenarios are characterized by their value of 
texp, ipo, and the stellar IMF. The IR luminosity (red curve) is characterized by its dust temperature. The dashed 
red line indicates the lower limit on ipo for the different SF scenarios (see Table 1). The filled circles on the two 
horizontal lines illustrate the connection between ipo and Td for the two cases shown in Figure 2c; (e)-(f) Comparison 
of modeled stellar and dust masses (hatched regions) and SED determined masses (diamonds) as a function of ipo 
(vertical ticks) and fn g (black lines). See text for details. 
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Fig. 3. — The evolution of the dust mass for a closed box model with a low-metallicity IMF and if ) o = 2 Mq yr _1 . 
The bold solid lines assume no grain destruction. The thin solid and dashed lines are models with different initial gas 
mass. The evolution of the gas mass determines the rate of grain destruction. Large gas masses imply low dust-to-gas 
mass ratios and lower rate of grain destruction [eqs. (5) and (6)]. The figure also shows that AGB stars do not 
contribute to the dust mass. 
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Fig. 4. — The wavelength dependence of k(A), derived from forcing the intrinsic spectra to reproduce the observed 
galaxy’s UVO spectrum, is compared to the mass absorption coefficient of 0.01 (jm radius astronomical silicate and 
amorphous carbon grains. The curves are labeled by the value of i/)o in units of Mq yr -1 . For a given optical depth 
and dust mass, k is proportional to the projected area of the disk. Plotted values assume an area of 1 kpc 2 . The 
derived mass absorption coefficient exhibits a similar rise at UV wavelengths as the small silicates grains. Young 
galaxies with a silicate dominated population of interstellar dust grains will therefore have a break in their spectrum 
around 0.15 ^m, which can be confused with the Lyman break. 
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Table 1. Viable star formation scenarios following successive model constraints 1,2 



Energy constraint 

Mass production constraint: Md(SEDfit) = Md(evol) 

SF scenario 

Labs = Lir 

M g (0) = 1 x 10 10 Mq 

M g { 0) = 5 x 10 10 Mq 

texp = 100 Myr Low-Z IMF 

ipo > 30 

ipo > 40 

ipo^40- 50 



M*> 7 x 10 s 

M*« (7-10) x 10 8 



M d < 6 x 10 7 

M d « (7 — 3) x 10 7 

Kroupa IMF 

ipo > 45 

no viable models 

no viable models 


texp = 300 Myr Low-Z IMF 

ipo > 4 

ipo » 5 - 35 

ipo « 5 - 8 



M (2 - 12) x 10 8 

M»^ (2 - 3) x 10 8 



M d « (70 - 5) x 10 6 

M d ss (7 — 2) x 10 7 

Kroupa IMF 

ipo > 10 

no viable models 

V>o >35 
M*> 1 x 10 10 
M d < 9 x 10 6 


1 Values of ip o are in units of Mg yr 1 2 . 

2 Values of the stellar and dust masses, M* and Md, respectively, and are in units of Mg. 


